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Abstract 

Direct and inverse Auger scattering are amongst the primary processes that mediate 
the thermalization of hot carriers in semiconductors. These two processes involve 
the annihilation or generation of an electron-hole pair by exchanging energy with a 
third carrier, which is either accelerated or decelerated. Inverse Auger scattering 
is generally suppressed, as the decelerated carriers must have excess energies higher 
than the band gap itself. In graphene, which is gapless, inverse Auger scattering 
is instead predicted to be dominant at the earliest time delays mm- Here, < 8 
femtosecond extreme-ultraviolet pulses are used to detect this imbalance, tracking 
both the number of excited electrons and their kinetic energy with time- and angle- 
resolved photoemission spectroscopy. Over a time window of approximately 25 fs 
after absorption of the pump pulse, we observe an increase in conduction band carrier 
density and a simultaneous decrease of the average carrier kinetic energy, revealing 
that relaxation is in fact dominated by inverse Auger scattering. Measurements of 
carrier scattering at extreme timescales by photoemission will serve as a guide to 
ultrafast control of electronic properties in solids for PetaHertz electronics [5j. 
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FIG. 1: Different Coulomb-interaction-mediated scattering processes in photo-excited graphene: 
impact ionization (green) and Auger heating (yellow). These processes can be identified experi¬ 
mentally by comparing the total number of electrons inside the conduction band, Nqb , with the 
average kinetic energy of electrons inside the conduction band, Ecb/Nqb- Occupied states are 
shown in red, empty states are shown in white. 


The dynamics of pre-thermal Dirac carriers in graphene are expected to host interest¬ 
ing and unconventional phenomena. Microscopic simulations ms predict that different 
scattering events mediated by the Coulomb interaction may contribute to the ultrafast re¬ 
distribution of the photo-excited electrons (see Fig. [Tj) . For example, these theories imply 
that inverse Auger scattering (also known as impact ionization) may dominate at early times. 

This can be understood as follows BE!. As shown in Fig. [lj for a non-equilibrium 
distribution with holes in the valence band at E = Ed — ump and electrons in the 

conduction band at E = Ed + \fajJpump (where Ed and hu pnmp are the Dirac point and 
pump photon energy, respectively) in undoped graphene, the recombination of electron-hole 
pairs required for Auger heating is strongly suppressed due to the lack of holes at the top 
of the valence band. On the other hand, the available phase space for the inverse process 
(impact ionization), where the excess energy of an electron high in the conduction band 
is used to generate secondary electron-hole pairs, is large. Therefore, impact ionization is 
believed to dominate over Auger heating for as long as it takes to establish a thermalized 
electronic distribution. The resulting carrier multiplication, for which the absorption of 
a single photon may generate multiple electron-hole pairs, has raised interest for possible 
applications in photovoltaic devices. 
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However, these theoretical arguments apply only for undoped graphene and low excitation 
fluences (few /iJ/cm 2 ). At high pump fluences, the balance between impact ionization and 
Auger heating is predicted to be reestablished within a few tens of femtoseconds, reducing 
the carrier multiplication factor considerably BE3. Furthermore, real graphene samples 
typically rest on a substrate resulting in a non-negligible doping of the graphene layer. 
The presence of either hole or electron doping will shift the relative importance of impact 
ionization and Auger heating, suppressing carrier multiplication. 

Previous optical experiments deduced transient electron distribution functions by com¬ 
paring differential transmission, reflectivity, or absorption data to model calculations. In 
this way, high-temporal-resolution experiments estimated an electronic thermalization time 
between 13 and 50 fs 0E]. Furthermore, indirect evidence for carrier multiplication with 
~200 fs pulses in the low fluence regime (< 30 /i.J/cm 2 ) was obtained by comparing the 
number of absorbed photons to the number of electron-hole pairs [4] with some indications 
already in 13]. 

Photoemission techniques, which directly measure electron numbers as a function of en¬ 
ergy and momentum, are ideally suited for a direct visualization of Auger scattering. Pre¬ 
vious time- and angle-resolved photoemission spectroscopy (tr-ARPES) experiments were 
performed at pump fluences on the order of mJ/cm 2 and thus short time scales for the 
initial thermalization. The temporal resolution of > 30 fs was then insufficient to resolve 
pre-thermal carrier distributions [MU- 

Here, we use tr-ARPES with < 10 fs pump and probe pulses to access charge carrier 
dynamics in photo-excited graphene. The output of a 1 kHz, 30 fs TitaniunuSapphire 
amplifier operating at 1.55 eV photon energy was spectrally broadened in a neon-filled 
hollow-core fiber and recompressed using a set of chirped mirrors, resulting in 8 fs pulses with 
a slightly blue-shifted spectrum. 400 /i,J of energy were converted into extreme ultraviolet 
(XUV) pulses by High-order Harmonics Generation (HHG) in Argon. A single harmonic at 
ha;p ro be = 30 eV was selected with a time-preserving grating monochromator [15| and used as 
a probe pulse for photoemission. The sample was excited with a fluence of 20 mJ/cm 2 with 
a pump pulse duration of ~ 10 fs. Both pump and probe pulses were s-polarized, with the 
electric held vector in the plane of the graphene sample perpendicular to the TK direction. 

Tr-ARPES experiments were performed on lightly hole-doped quasi-freestanding epitaxial 
graphene monolayers with the Dirac point 200 meV above the equilibrium chemical potential 
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FIG. 2: tr-ARPES measurements. A Snapshots for different pump-probe time delays. B Cor¬ 
responding pump-induced changes of the photocurrent. The data in panel A and B has been 
smoothed. C Electronic distribution functions (red line) together with Fermi-Dirac fits (FD fit, 
black dashed line) revealing the presence of non-thermalized carriers (NTC, blue). The y-axis scale 
in C is linear. The measurements presented in this figure were carried out using the high-energy- 
resolution grating G300 (details see text). 


[lj]j. These measurements were taken with two different monochromator gratings, optimizing 
either time or energy resolution [[15] . 

In Fig. [2]A we show tr-ARPES snapshots of graphene’s linear 7r-bands along the TK 
direction for selected pump-probe time delays across the rising edge of the pump-probe 
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signal. Due to the s-polarization of the probe pulse ra and to photoelectron interference 
effects PI] only the right-hand branch of the Dirac cone was visible in this geometry. These 
snapshots were recorded with a temporal resolution of 14 fs and an energy resolution of 500 
meV (details see Supplementary Information) achieved by using a monochromator grating 
with 300 grooves per millimeter (G300) pSJ- The optical matrix element describing the 
absorption of the pump photon is anisotropic with nodes along the direction of the pump 
polarization (in this case perpendicular to the TK direction) and maxima in the direction 
perpendicular to the pump polarization (in this case TK) jTTTl f?T] . We probed the response 
of the electronic structure along the TK direction where the effect of the pump pulse is the 
strongest. 

The pump-induced changes of the photocurrent are plotted in Fig. [2]B. A loss (gain) of 
electrons below (above) the equilibrium chemical potential, which is used here as zero-energy 
reference, is observed. 

From the snapshots in Fig. [2]A, transient electron distribution functions were extracted 
and compared to Fermi-Dirac (FD) distributions in Fig. [2p. At negative delays, immediately 
before arrival of the pump pulse (-10 fs in Fig. [2p), the distribution follows a FD distribution, 
indicating the presence of a completely thermalized electron gas. Near zero time delay, a 
shoulder develops in the conduction band above Ed = 200 meV. This shoulder cannot be 
fitted with a FD distribution. The residual weight between the experimental distribution 
function (red line in Fig. |2p) and the FD fit (dashed black line in Fig. [2p) is taken as 
proportional to the number of non-thermal carriers (NTCs, blue shaded area in Fig. IP)- 
We find that the number of NTCs keeps increasing until the peak of the pump-probe signal 
is reached at a time delay of ~20 fs. These NTCs can be understood as a precursor of the 
inverted carrier population observed in previous studies with longer pump pulses at slightly 
smaller photon energies |8j, [13j [22] . 

Higher-temporal-resolution measurements were performed with a second monochromator 
grating with 60 grooves per millimeter (G60), delivering in a temporal resolution of 8 fs 
and an energy resolution of 800 meV (details see Supplementary Information). The num¬ 
ber of carriers in the conduction band at E > Ed, Nqb , and their average kinetic energy, 
Ecb/Ncbi were determined directly from the raw data (details see Methods and Supple¬ 
mentary Information) and displayed in Fig. [3j For about 25 fs around zero pump-probe 
time delay, the average kinetic energy Ecb/^cb was observed to decrease while the number 


5 


T- 1.0 



to 

< 


h 0.2 


FIG. 3: Direct evidence for impact ionization. Comparison between the temporal evolution of the 
total number of carriers inside the conduction band ( Ncb , light red) and the temporal evolution 
of their average kinetic energy ( Ecb/Ncb , light blue). Around zero time delay Ecb/Ncb already 
decreases while Ncb keeps increasing, indicating impact ionization. Dark red and blue lines are 
fits to the data that serve as guides to the eye. The fitting function consists of an error function 
to describe the rising edge plus a single (Ncb) or double exponential decay (Ecb/Ncb)- The 
measurements presented in this figure were carried out using the low-energy-resolution grating 
G60 (details see text). Error bars represent the standard deviation. 


of carriers Ncb kept increasing, indicating impact ionization. 

From these measurements the following scenario can be envisaged. Impact ionization 
is the primary scattering mechanism during the first ~25 fs, accumulating carriers at the 
bottom of the conduction band and establishing a precursor of the population inversion 
observed previously mmm- This state then likely decays through Auger heating and 
electron-phonon scattering within ~100 fs, re-establishing a single Fermi Dirac distribution 
[8j. We note that the temporal evolution of the electronic temperature and the decay time of 
the non-thermal carriers found in our experiment (see Supplementary Information) further 
substantiate this interpretation. 

In summary, we have used tr-ARPES with < 10 fs pulses to identify the primary scat¬ 
tering events that result in the rapid thermalization of photo-excited electron-hole pairs in 
graphene. By comparing the number of electrons inside the conduction band with their 
average kinetic energy we find that impact ionization is the predominant scattering chan¬ 
nel within the first ~25 fs also in doped graphene at high fluence. Whether the observed 
carrier multiplication can be exploited for solar cell applications remains questionable, as 
the absence of a band gap in graphene makes charge separation difficult. Nevertheless, the 
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ultrafast dynamics of photo-excited electron-hole pairs observed here will be of use for the 
design of new electronic and optoelectronic devices operating at Pffz rates [3] based on 
graphene and other materials. 


METHODS: 

Samples. Quasi-freestanding epitaxial graphene samples on silicon carbide were grown 
as previously described in [16] . Prior to graphene growth the silicon carbide substrates were 
etched in hydrogen atmosphere to remove scratches from mechanical polishing. In a second 
step the substrates were annealed in argon atmosphere resulting in the growth of one car¬ 
bon monolayer on the silicon-terminated face of the substrate. This carbon monolayer was 
subsequently decoupled from the substrate by hydrogen intercalation and characterized by 
static ARPES measurements. The samples were transported to the Artemis user facility un¬ 
der ambient conditions, reinserted into ultra-high vacuum, and cleaned by a mild annealing, 
recovering the original band structure. 

Tr-ARPES. Tr-ARPES experiments were performed using the Materials Science end 
station at the Artemis user facility at the Rutherford Appleton Laboratory in Harwell, 
United Kingdom. The setup consists of a TitaniuimSapphire amplifier (hu = 1.55 eV) 
operating at 1 kHz with 30 fs pulse duration. 1 mJ of energy was compressed down to 
8 fs with a hollow-core fiber filled with 1.4 bar of neon operated with a pressure gradient 
and followed by 10 chirped mirrors. Pulse compression resulted in a slight blue shift of the 
spectrum. 400 /i J of energy were used for high harmonics generation (HHG) in argon, 20 
mJ/cm 2 were used to excite the sample. Due to the longer path in air two additional chirped 
mirrors were added to the pump path resulting in a pump pulse duration of 10 fs. All pulse 
durations were measured using second-harmonic frequency-resolved optical gating (FROG). 
Out of the broad HHG spectrum one particular energy, hoj pio b e = 30 eV, is selected by a 
time-preserving grating monochromator [15] and used as a probe pulse for photoemission. 
For the present investigation two different gratings with different dispersion (G300 with 300 
grooves per mm and G60 with 60 grooves per mm) have been used for wavelength selection 

m- 

Data analysis. The data in Fig. [2]A and B has been smoothed. To obtain the electron 
distribution functions in Fig. If we took lineouts at constant angle (energy distribution 


7 


curves) from the raw data, determined their integrated intensity within the full width at 
half maximum (FWHM) of the peak and attributed the resulting number to the energy of 
the peak position. The number of carriers in the conduction band in Fig. 3 was obtained 
by integrating the electron distributions for E > Ed. The kinetic energy of carriers in the 
conduction band in Fig. [3] was obtained by integrating the electron distributions multiplied 
by the energy measured with respect to the Dirac point for E > Ed- 
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SUPPLEMENTARY INFORMATION: 


Sample characterization 


After growth and hydrogen intercalation the graphene samples were characterized using 
high-resolution angle-resolved photoemission spectroscopy (ARPES). In Fig. [4] we show the 
electronic structure along a cut through the K point perpendicular to the TK direction. The 
sample is lightly hole-doped with the Dirac point ~200 meV above the equilibrium chemical 
potential. 



FIG. 4: High-resolution ARPES spectrum of hydrogen-intercalated monolayer graphene on silicon 
carbide measured along a cut through the K point perpendicular to the TK direction. 


Pulse compression 

Pulse compression from 30 fs down to 8 fs was achieved by broadening the output spec¬ 
trum of a Titanium:Sapphire amplifier in a 1.4 bar neon-filled hollow-core fiber operated 
with a pressure gradient and subsequent recompression using 10 chirped mirrors. The pulse 
duration was measured with second-harmonic frequency-resolved optical gating (FROG). 
For that purpose the pulse was split by a 50/50 beam splitter, and the two beams were over¬ 
lapped again on a BBO crystal. The generated second harmonic signal was analyzed with a 
spectrometer. Figure |5]A shows the two-dimensional FROG trace acquired as a function of 
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FIG. 5: Second-harmonic FROG reconstruction of the compressed pulses. A Second-harmonic 
FROG trace acquired as a function of delay. B Reconstructed second-harmonic FROG trace. C 
Reconstructed pulse intensity profile. The retrieved pulse duration is 8.3 fs. 


delay between the two pulses. The reconstructed FROG trace and pulse intensity profile are 
reported in panels B and C of the same figure, respectively. The retrieved pulse duration 
is 8.3 fs. The delay-integrated spectrum is centered around 335 nm (3.70 eV), indicating a 
central wavelength of the compressed pulse of 670 nm (1.85 eV). 

This beam was split with an 80/20 beam splitter and used for high harmonics generation 
(HHG) in argon (400 /i.J) and to excite the sample (3.5 //J), respectively. Due to the longer 
path in air, two more chirped mirrors were added to the pump path, resulting in a pump 
pulse duration of 10 fs. 

An estimate for the temporal resolution of the tr-ARPES experiment can be obtained 
from the rise time of the pump-probe signal. For that purpose, the rising edge of the number 
of non-thermal carriers (see below) has been fitted with an error function, erf((f—f 0 )/rt). The 
rise time, rt, is 14 ± 2 fs for grating G300 (300 grooves per mm) and 8 ± 1 fs for grating G60 
(60 grooves per mm). This rise time is related to the full width at half maximum (FWHM) 
of the cross correlation between pump and probe pulses via FWHM = 1.665 x rt. The 
energy resolution has been obtained by fitting electron distribution functions (see below) at 
negative delay with a Fermi-Dirac distribution convoluted by a Gaussian. The temperature 
was fixed to 300 K (the temperature of the sample). The FWHM of the Gaussian was 
obtained as a fitting parameter, giving FWHM = 500 meV for G300 and FWHM = 800 
meV for G60. 
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Extraction of electronic distribution functions 


The energy resolution of the tr-ARPES measurements was restricted to 500 meV using 
grating G300 in Fig. [2] and 800 meV using grating G60 in Fig. 3 of the main manuscript. 
This means that the established standard methods to extract electron distribution functions 
(simple momentum integration |8j or MDC method described in [23] ) fail, meaning that 
the equilibrium distribution function determined that way does not follow a Fermi-Dirac 
distribution. 

Our method is based on the assumption that the comparatively low energy resolution 
in our experiment broadens the data along the energy axis, while the angular resolution is 
unaffected. Thus, we took lineouts at constant angle (energy distribution curves, EDCs), 
determined their integrated intensity within the FWHM of the peak and attributed the 
resulting number to the energy of the peak position. This procedure is illustrated in Fig. [6} 
The resulting distribution functions are shown in Fig. [2p of the main manuscript. 

From similar distribution functions recorded using G60 we calculated the number of 
carriers inside the conduction band, N C b, and their average kinetic energy, E C b/Ncb , as 
described in the following. In Fig. [6p the electronic distribution is plotted as a function of 
energy and pump-probe delay. Due to the light hole-doping the bottom of the conduction 
band, Ed, lies ~200 meV above the equilibrium chemical potential, /i e . Ncb is determined 
by summing up the intensity between E = Ed and E = 2.6 eV for a given delay. Similarly, 
Ecb is determined by summing up the intensity multiplied by the energy referenced with 
respect to the Dirac point (Ed = 0 eV) between E = Ed and E = 2.6 eV for a given delay. 

Results of Fermi-Dirac fits 

Figure [7] summarizes the results of the Fermi-Dirac fits from Fig. IP in the main 
manuscript. The figure shows the temporal evolution of the electronic temperature (T e , 
light red), and of the number of non-thermal carriers (NTC, light blue), where NTC corre¬ 
sponds to the area of the blue shaded region in Fig. [2p of the main manuscript. Dark red 
and blue lines represent fits to the data including an error function to describe the rising 
edge and a single exponential decay. 

The temperature peaks 30 fs after the number of non-thermal carriers reaches its max- 
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FIG. 6: A Electron distributions in B were extracted from the raw data by summing up the 
intensity along the white arrow and attributing it to the energy indicated by the black arrow. 
Blue lines indicate the peak position (continuous) ± full width at half maximum (dashed). B 
Electron distribution function as a function of energy. The drop at negative energy is caused by 
photoemission matrix element effects. C Electron distribution functions from B for different pump- 
probe delays as a two-dimensional color plot. The position of the chemical potential fj, e and the 
Dirac point Ed are indicated by white and black lines, respectively. 


imum. The decay times of the temperature and of the number of non-thermal carriers are 
412 ± 32 fs and 113 ± 6 fs, respectively. In graphene, the elevated electronic temperature is 
known to cool down by optical and acoustic phonon emission [2D - l2fi] . The life time of the 
non-thermal carriers is similar to the one found for the population inversion in [8, 13] (22J, 
indicating that complete thermalization is hindered by the relaxation bottleneck imposed 
by the Dirac point. 


tr-ARPES snapshots for grating G60 

The tr-ARPES data recorded with grating G60 that has been used to generate Fig. [3] in 
the main manuscript is displayed in Fig. [8j Similar to the data presented in Fig. [2]A and B 
the data in Fig. [3] has been smoothed. Data displayed in Fig. [2p and Fig. [3] in the main 
manuscript and Fig. [7] of the Supplementary Information was extracted directly from the 
un-smoothed raw data as explained above. 
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FIG. 7: Results of Fermi-Dirac fit of the data shown in Fig. [2p of the main article. Electronic 
temperature (T e , light red) and number of non-thermalized carriers (NTC, light blue) as a function 
of pump-probe delay. Zero time delay is defined as the middle of the rising edge of the number 
of non-thermalized carriers. Dark red and blue lines are fits to the data. Error bars represent the 
standard deviation. 
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FIG. 8: tr-ARPES snapshots recorded with grating G60. 
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